OTHER KUIPER BELTS 



M. Jura 

Department of Physics and Astronomy, University of California, Los Angeles CA 

90095-1562; jura@clotho.astro.ucla.edu 

ABSTRACT 

When a main sequence star evolves into a red giant and its Kuiper Belt 
Object's (KBO's) reach a temperature of ~170 K, the dust released during the 
rapid ice-sublimation of these cometary bodies may lead to a detectable infrared 
excess at 25 /im, depending upon the mass of the KBO's. Analysis of IRAS data 
for 66 first ascent red giants with 200 Lq < L < 300 L© within 150 pc of the Sun 
provides an upper limit to the mass in KBO's at 45 AU orbital radius that is 
usually less than ~0.1 M^. With improved infrared data, we may detect systems 
of KBO's around first ascent red giants that are analogs to our Solar System's 
KBO's. 

Subject headings: circumstellar matter - comets - stars, red giants 



1. INTRODUCTION 

The formation and evolution of large solids such as comets, asteroids and planets in 
astrophysical environments is of great interest. Here, we investigate whether other stars 
possess Kuiper Belt Objects (KBO's) similar to those found in the Solar System. 

We assume that the hypothetical KBO's around other stars resemble those in our Solar 
System and have a composition of ice and dust similar to that of comets (see Luu & Je- 
witt 2002). During the main sequence phase of the star's evolution, the KBO's quiescently 
remain in stable orbits. However, when the star becomes a red giant, the KBO's become 
sufficiently hot that ice sublimates and previously embedded dust particles are ejected into 
the surroundings. In §2, we sketch our argument, which is presented in full detail in §3, 
that sufficient dust may be released during this red giant phase that the star can display a 
detectable infrared excess. 

Stern, ShuU & Brandt (1990) and Ford & Neufeld (2001) computed the fate of ice 
sublimated from comets that have orbital radii greater than 100 AU when the host star 



- 2 - 



evolves onto the Asymptotic Giant Branch and attains a luminosity near 10'' Lq. Now that 
the Kuiper Belt in the Solar System with most objects lying near 45 AU from the Sun is 
becoming better understood (Luu & Jewitt 2002), it is possible to imagine the response of 
such a system to a star's first ascent up the red giant branch when its luminosity exceeds 
100 L0. However, although progress is being made, there are still significant uncertainties 
about the Kuiper Belt in the Solar System. Recent estimates of its mass range from 0.01 
Me (Bernstein et al. 2003) to 0.1 (Luu & Jewitt 2002). Therefore, although we describe 
a procedure to investigate Kuiper Belt-like systems around other stars, our own outer Solar 
System is so poorly understood that an exact comparison is not yet possible. 

Two current descriptions of the KBO's differ mainly but not exclusively in estimates of 
the numbers of objects with radii smaller than ~200 km. With the expectation that short 
period comets originate in the Kuiper Belt, Luu & Jewitt (2002) propose a rapid rise in the 
numbers of such smaller KBO's. In contrast to this expectation, Bernstein et al. (2003) used 
the Hubble Space Telescope to search for KBO's as faint as 28.3 mag, corresponding to radii 
of ~15 km, and found many fewer objects than expected from the model given by Luu & 
Jewitt (2002). However, with this observational result of Bernstein et al. (2002), the origin 
of the short period comets becomes a mystery. Perhaps there is a bimodal distribution of 
KBO's. In this paper, given the uncertainties in current knowledge of the outer Solar System, 
we consider models both with and without large numbers of KBO's smaller than 200 km in 
radius. 

Most first ascent red giants do not show any evidence for circumstellar dust. The 
minority of first ascent red giants with infrared excesses (Judge, Jordan & Rowan- Robinson 
1987, Jura 1990, 1999, Kim, Zuckerman & Silverstone 2001, Smith 1998, Zuckerman, Kim & 
Liu 1995, Plets et al. 1997) do not seem to have sublimating KBO's, since the characteristic 
inferred dust temperature is less than 70 K instead of 170 K, as we predict in §3. Here, we 
focus on the possibility of detecting relatively warm dust from disintegrating KBO's. 

2. OVERVIEW 

Because the details of our model are somewhat complicated, we first present a schematic 
overview of our analysis, while the complete description is presented in §3. Using our own 
Kuiper Belt as a prototype, we assume that the KBO's orbit the host star near 45 AU. At this 
orbital separation, we find that most of the KBO material survives in the solid phase as long 
as the temperature is less than ~150 K which occurs as long as the star is less luminous than 
~170 Lq. However, once the star's luminosity exceeds this value, the KBO's become warm 
enough that relatively rapid ice sublimation occurs. In the time interval, approximately 
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characterized by the parameter tjioB, during which the star increases in luminosity from 
~170 L0 to ~300 Lq, most of the mass of the KBO's is lost into their surroundings. As 
the KBO's are sublimating, some dust particles are released that continue to orbit the star. 
If sufficient mass is released from the KBO's, these orbiting dust grains can absorb enough 
light to produce a detectable infrared excess. 

If MKBoi^) denotes the total initial mass of the KBO's and forbu denotes the fraction of 
this mass which is ejected as dust particles that are large enough to remain in orbit around 
the star and are not blown out of the system by radiation pressure, then the total amount 
of mass that can orbit the star is forbuMKBoi^)- We assume that the dust particles are 
spheres of radius b and density Pdust, and they orbit the star of luminosity, L^, at distance, 
D. While this assumption is relaxed in §3.7, here, for simplicity, we assume that D does not 
change during a dust particle's evolution. Also, while other factors are considered in §3.7, 
here, we assume that the lifetime of the orbiting dust particles is controlled simply by the 
Poynting-Robertson time, tpR. Following Burns, Lamy & Soter (1979), we write: 



Since the Poynting Robertson lifetime is typically shorter than the the characteristic 
evolutionary time on the red giant branch, at any given moment during the star's evolution- 
ary phase when grains are being released, the fraction of orbiting material contributing to the 
infrared excess is tpn/tnGB- If the grains have opacity, x (cm^ g^^), and if the circumstellar 
dust cloud is optically thin, then the total luminosity of the circumstellar dust, L excess is 
scaled from the fractional dust coverage of an imaginary sphere of surface area AtiD"^ by the 
relationship: 

. _ ( Xf orbit MkBo{^) \ ( tpR\ 

~ I ; \j^b) * 

We use a simple model of the grain opacity where the geometric cross section (7r6^) equals 
the total absorption cross section. In this case, since each dust grain has mass (47rp(iust^'^/3), 
we may write that: 

^PdustO 

Combining equations (l)-(3), we find: 

J _ f orbit MkBo{0) ... 

^excess ~ i , 

^tRGB 

Expression (4) captures the essence of our model. The infrared luminosity simply scales as 
the rate of production of dust grains which varies as MKBo{0)/tRGB- If f orbit and Mkbo{0) 
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are sufficiently large, then the infrared excess may be detectable. Our estimate of Lexcess, 
which is computed more exactly in §3, is independent of the size of the KBO's as long as 
these objects are fully destroyed by sublimation. 

In §3, we find characteristic values of forbit ~ 0.25 and tnoB ~ 14 Myr. Therefore, if 
Mkbo{Q) ~ 0.1 M®, then from equation (4), Lexcess ~ 0.02 Lq. Because the KBO's mostly 
sublimate near 170 K, the spectrum of this excess luminosity peaks at wavelengths near 25 
/im. In §4, we show that a 170 K black body with a luminosity of 0.02 Lq is just at the 
threshold of having been detected with IRAS data for red giants with total luminosities less 
than 300 L^. Consequently, the absence of a measurable infrared excess at 25 fim. around 
these red giants can be used to place an upper limit to the mass of KBO's closer than ~50 
AU of ~0.1 M®. In §5, we discuss our results and in §6, we present our conclusions. 



3. DETAILED MODEL 

We now elaborate about the schematic model described in §2. To compute the fate of 
the KBO's when a star is evolving as a red giant, we write for the total mass loss rate from 
KBO's, Mkbo, that: 

MkBO = CTKBoiT) AKBoit) (5) 

In equation (5), &kbo denotes the mass loss rate per unit surface area of the KBO's as a 
function of temperature, T, while AxBoit) denotes the total surface area of these objects 
as a function of time, t. In §3.3 and §3.4, we describe calculations for &kbo and Akbo as 
functions of time or temperature, and we show in §3.7 how Mrbo can be used to estimate 
the infrared excess. 



3.1. Initial Mass and Total Area of the KBO's 

Here, we consider two size distributions: systems with initially many small (radius less 
than 200 km) KBO's, and Kuiper belts with initially no KBO's smaller than 200 km in 
radius. 
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3.1.1. Model with initially many small KBO's 

Following Luu & Jewitt (2002), we assume the KBO's are icy spheres of radius, a, and 
follow a broken power law for the size distribution, n(a) da. For a > abreak, 

n{a) da — no ^ ^^^e^fc j 

while for a < abreak, 

n{a) da — no ^ ^^^"■'^ ^ ^^Jq 

From Luu & Jewitt (2002), we take af,reak — 1 km and a maximum size, a^ax, of about 
1000 km. The distribution in equation (7) for the KBO's smaller than 1 km is derived from 

theoretical models (Kenyon 2002) and is consistent with upper limits to the contribution 
of KBO's to the zodiacal light and to distortions of the microwave background (Kenyon & 
Windhorst 2001, Teplitz et al. 1999). 

Before a star becomes a red giant, the initial mass of the larger KBO's, Miargei^), is: 

Miarge{0) = / Tlo j I 1 da = — Pi^fiO ^Iq Ofcreaifc ( ) (8) 

^ ^break ^ \ / \ ^break / 

where Pkbo denotes the density of the KBO's. Following Greenberg (1998) and Kenyon 
(2002), we assume that the KBO's are composed mostly of ice with some additional refractory 
dust and that pkbo = 1-5 g cm~^. From the measured size distribution of KBO's (Luu & 
Jewitt 2002), we adopt Miargd^) = 0.1 as a first approximation to the mass of KBO's. 
Luu & Jewitt (2002) use Pkbo — 1-0 g cm~^ and a shghtly smaller total mass of 0.08 M^. 

The total initial surface area of the large KBO's, Aiargd^), is 
AargeiO) = Uo (^) ' (4 TTa^) da ^ 47rnoaL„. «^ ^ (9) 

J abreak ^ pKBO abreak In 

For the parameters we adopted above, ^iarge(O) = 1-7 x 10^^ cm^. 

We can also compute the total mass contained in the KBO's with a smaller than 1 km. 
If dmin denotes the smallest KBO, then: 

/n^ T'""' fabreak^-^ (^T^PKBOa^\ , SvT ^ 

MsmaiiK^) = j ^^\^) V 3 ) ~ Y^^so no a^^g^fc (10) 

As long as amin « abreak, then evaluation of the integral in equation (10) shows that the 
mass of the small KBO's is essentially independent of a^j„. From equations (8) and (10), 
we find that: 
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M™,,(0) ^ (11) 



In ( ) 



Therefore, with our adopted values of amax ^-nd Qbreak, we find that Msmaiii^) ~ 0.3 Miargei^)', 
most of the mass of the KBO's is contained in the larger objects. 



3.1.2. Model with initially only large KBO's 

As a second model, we adopt a simplified version of the results of Bernstein et al. (2003). 
In particular, we assume that the size distribution of KBO's follows equation (6), but there 
are no KBO's smaller than Ubreak with abreak — 200 km instead of 1 km. We denote the initial 
mass of the KBO's in this system as Mbig{0). If the numbers of objects larger than 200 km 
in radius is the same as in the model described in §3.1.1, then Miyig{0) ~ 0.02 M^. Although 
this model is representative of the results of Bernstein et al. (2003), it does not agree exactly 
with their description of the Kuiper Belt which involves either a "rolling power law" fit or 
a "double power law" fit to the size distribution of the KBO's. One source of uncertainty 
is that even for KBO's larger than 200 km in radius, there is some difference between the 
descriptions of Luu & Jewitt (2002) and Bernstein et al. (2003) which we cannot resolve in 
this paper. In any case, by evaluating the integral in equation (9) with the new lower limit 
to the size of 200 km, we find that the initial area of the KBO's in this model, Aing{0), is 8.7 
X 10^^ cm^. 



3.2. Red Giant Branch Evolution 

The fate of a KBO is controlled by the host star's evolution. For simphcity, we param- 
eterize the detailed calculations by Girardi et al. (2000) for a star's luminosity, L*, on the 
red giant branch as: 

L,(i) RiLoe*/*«°« (12) 

For stars of 1 Mq, we fit equation (12) with tRCB — 1-4 x 10^ yr, Lq = 87 Lq and a starting 
time {t — 0) of 1.220 x 10^° year after the star's birth. For stars with 1.5 Mq, we employ 
the same value of tROB and adopt a starting time of 2.7742 x 10^ yr and Lq = 91 Lq. In 
Figure 1 we show a comparison of the prediction from equation (12) using tRCB = 1-4 x 10^ 
yr with the detailed calculations by Girardi et al. (2000) for stars of 1 Mq and 1.5 Mq. For 
luminosities less than 1000 Lq, the agreement is mostly better than 20%. Although equation 
(12) is an imperfect match to the detailed calculations, it docs illustrate that the growth of 
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luminosity with time on the red giant branch for a star of ~ 150 L© is very approximately 
exponential with an e-folding time of ~ 1.4 x 10^ yr. 



3.3. Ice Sublimation 

With the description given in §3.2 of the star's evolution, we can estimate the rate of 
ice sublimation from an individual KBO. We start with the expression: 

^ g = -47ra a KBO (13) 

Therefore, at time tf, the total decrease in the radius of a KBO, Aa, is given by the following 
solution to expression (13): 

A„(t,) = /" WIW) (14) 
J Pkbo 

As in equation (5), T denotes the grain temperature. Note that Aa is independent of a. 

To estimate axso we follow Ford & Neufeld (2001) who denote this quantity as rh. 
Converting from the units they employ and using their expression for the vapor pressure of 
ice, we write that 

ctkbo = aoT-V^e-^-"/^ (15) 

where (Tq = 3.8 x 10^ g cm~^ s~^ K^^'^ and Tguu — 5530 K. In expression (15), the chemical 
binding energy of H2O molecules onto the solid ice is ksTsuhi where ks is Boltzmann's 
constant. 

To compute Aa from equation (14), we need to find the temperature of a KBO as a 
function of time. If we ignore the difference between the day an night side of the KBO, 
assume unit emissivity and an albedo of 0, close to 0.07 which may be a typical for a KBO's 
albedo (Luu & Jewitt 2002), then for a KBO at distance, -Dinit, from the star: 

Tit) = V^' (16) 

where asB is the Stephan-Boltzmann constant. In equation (16), we ignore cooling by subli- 
mation. At 170 K, the characteristic temperature when most of the ice has been sublimated, 
it can be shown from the results in Ford & Neufeld (2001) that cooling by ice sublimation 
is less than 7% of the cooling by radiation. 

From equations (12) and (16), we write that 

djT 

dt = Urgb^ (17) 
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Using equations (15) and (17), we re-write equation (14) to find: 

Aa{Tf) = 4^^^ p T-^/^ ^-Tsu,UT (18) 
Pkbo J 

where Tj is the KBO's temperature at time tf. With the substitution that v? = Tguu/T, the 
integral in equation (18) can be re-expressed as a complementary error function to give: 



oo 

2 



Aa(T,) = 8 ^^^^j / e-«^ du (19) 

PKBOJ-subl-' 

When Uf > 1 which is equivalent to our case where Tsuu » T, then (see, for example, 
Abramowitz & Stegun 1965) 



/ 



„2 



le-"/ 



e"" du^ (20) 



Therefore: 



Aa{Tf) ^ 4 ^ e-^-«/^/ (21) 

Pkbo Tguu 

Equation (21) describes the shrinkage of KBO's as a function of their temperature. Numer- 
ically, for Tf — 150 K, 160 K, and 171 K, then Aa = 1 km, 10 km, and 100 km, respectively. 
These temperatures are achieved by KBO's at 45 AU when the red giant's luminosity is 170 
L0, 220 L0, and 290 Lq, respectively. Complete destruction of the KBO's is achieved when 
Aa = 1000 km which occurs for KBO's at 45 AU, when the star's luminosity is ~400 Lq. In 
the model with initially many small KBO's at 45 AU orbital radius described in §3.1.1, by 
the time the temperature attains ~170 K, all the KBO's smaller than 100 km are destroyed 
and somewhat more than half of the total initial mass of KBO's has been subhmated. There- 
fore, we adopt 170 K as a characteristic reference temperature for production of dust from 
KBO's. 



3.4. Evolution of the Total Area of the KBO's 

Using our calculation for the evolution of the size of a single KBO, we now consider the 
evolution of the surface area of the entire ensemble to use in equation (5). Each KBO of 
initial radius a evolves into a KBO with radius a' where 



a' — a — Aa(i) 



(22) 
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3.4-1- Model with initially many small KBO's 
For the large KBO's, we use equation (6) to find that their size distribution, n(a') da', 

n{a')da' = no (^r?p^) (23) 

To evaluate the total area of the KBO's with time, there are two regimes to consider. If Aa 
< flfcreofe, we write that 

(47ra'^)noU^5^ da' (24) 

a,reak-Aa \a' + Aa J 

Using equation (9), we evaluate expression (24) to find in the limit where a^ax >> O'break 
that: ^ 

where 5 — Aa/ahreak and Aiarge{^)- Prom equation (25), we see that the total surface area 
of the large KBO's remains nearly constant until Aa approaches abreak- Note that equation 
(25) is valid only when S < 1. 

When Aa > a^reak we must use a different lower bound in the integral and we write that 

(47ra'^)no(^^;^^j da' (26) 
We evaluate equation (26) to find in the limit where amax » Aa that: 

Alargeit) ~ Aiargei^)-^^^ = (2''') 

Equation (27) is valid when 6 > 1. In this regime, the total area of the large KBO's decreases 
as Aa increases with time. As Aa approaches amax, equation (27) overestimates Aiarge{t)- 
However, as long as Aa < 0.5 amax, then we can find from equation (26) that equation (27) 
is accurate to better than a factor of two. 

In the above paragraph, we described the evolution with time of the total surface area 
of the KBO's initially larger than abreak- In this paragraph, we describe the evolution of the 
surface area of the smaller KBO's. As long as Aa > a^j„, we can write from equation (7) 
that: 

(4 7ra'^)nof^;^^J da' (28) 
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Expression (28) can be conveniently evaluated to give: 

Asmaiiit) = (2 (1 - 5^1-) - 1(1 - 5^1-) + |(1 - 5'/^)) (29) 

2pKBo{abreak^a)'-/^ \ 3 5 / 

Note that Agmaii approaches zero as Aa approaches af,reak, or, equivalently, as S approaches 
unity. 

3.4-2. Model with initially only large KBO's 

To calculate Ai,ig{t) we simply substitute for Aiaj.gf.{t) in the equations in §3.4.1 with 
O'break = 200 km instead of 1 km. When Aa < af,reak, we use equation (25), and when Aa > 
O'break-i wc usc equation (27). 

3.5. Total KBO Mass Loss Rate 

Armed with the KBO temperature and total surface area as functions of time, we can 
evaluate equation (5) to find the total rate at which KBO's loss mass into their surroundings. 
It is convenient to present Mkbo as a function of the host star's luminosity, since that is 
a measurable quantity. The star's evolution is given by equation (12), and the sublimation 
rate from the KBO's is given by equation (15). 



3.5.1. Model with initially many small KBO's 

For systems with many small KBO's, we consider two regimes for Mkbo as a function 
of the star's luminosity. First, when Aa > abreak, ^iarge(^) from equation (27) gives the total 
surface area of the KBO's. In this case, Mkbo reaches a saturation value and is nearly 
independent of L^. At saturation, we find from equations (5), (9), (15), (21) and (27) that 

4 In I .9max_ \ j \ I J \ tRGB J 

y O'break J / 

Equation (30) shows that when when Aa > abreak, the sublimation rate from the KBO's 
depends mainly upon the initial mass of the large KBO's and the time scale of the star's 
evolution on the red giant branch. In this phase, the rapid rise in the sublimation rate with 
temperature given in equation (15) is roughly balanced by the decrease in the total area of 



Mkbo{T) 
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KBO's given in equation (27). When T = 170 K, we find from equation (30) that Mrbo ~ 

1.2 Miarge{0)/tRGB: Or 1.6 X lO^^ g g-l. 

A second regime to consider is when Aa < a^^ak- In this case, the total surface area of 
the KBO's is given by the sum of Aiarge{t) and Asmaii{t) from equations (25) and (29). In 
this regime, Mrbo starts low and rises rapidly with the star's luminosity until saturation is 
achieved. 

Graphical results for Mrbo vs. luminosity and presented in Figure 2 for a model with 
Dinit — 45 AU and Mi^rgei^) = 0.1 M^. For comparison, we also show the results for models 
with Dinit of 40 AU or 50 AU. In all three cases, we see that Mrbo approaches the saturation 
rate predicted by equation (30) for sufficiently large value of the star's luminosity. 



3.5.2. Model with initially only large KBO's 

If initially there are only large KBO's, then for much of the time during the star's 
evolution the total surface area of the KBO's is approximately constant. As long as Aa < 
200 km, we may write from equations (5). (15) and (25) that: 

Mrbo{T) = A,,g{0) (^1-6+ doT-i/^e-^-«/^ (31) 

Graphical results for Mrbo vs. luminosity from equations (30) and (31) are shown in Figure 
2 for the case where Dinn = 45 AU. We see that for stellar luminosities less than ~ 300 Lq, 
the dust production rate is relatively small in this regime where initially there are only large 
KBO's, but as the luminosity increases, the saturation value of Mrbo is attained. 



3.6. Dust Production 

In §3.5, we estimated the rate at which ice subhmates from the KBO's; here, we estimate 
the fractional mass of the KBO's released in dust that continues to orbit the star, f orbit- It 
is this orbiting dust which controls the infrared excess of the system. 

It is plausible that half of the total mass of the KBO's is dust (Luu & Jewitt 2002), and 
forbit could be as large as 0.5. However, we expect that f orbit is less than 0.5 since dust grains 
that are smaller than a critical radius, bcrit, do not orbit the star of mass, M*, and instead 
are driven out the system by radiation pressure. The critical size for being gravitationally 
bound is (Artymowicz 1988): 

b / 

""^ IGTTpdustGM^C 
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If the dust is composed of silicates with a density, Pdust, of 3 g cm ^, then for a 1 Mq star 
of 300 L0, the maximum luminosity that we consider, bcru — 57 /im. 

We assume that the size distribution of dust ejected from KBO's, Udustib) dh, resembles 
that of comets (Manner et al. 1992, Marker et al. 2002, Li & Greenberg 1998) which we 
approximate as: 

ridust = nc dh (33) 

We use hjnin and h^ax to denote the minimum and maximum size, respectively, of the dust 
from the KBO's while nc is a normalization constant which does not matter for our purposes 
here. We may write that the total mass in the dust, M^otai is 



total 



r-('-^)ncl-'dl (34) 

J brain V / 



Similarly, if Morut denotes the mass of the dust which is large enough to remain in orbit 
after it is produced, then 



MorUt = / I 7, ) T^co dh (35) 

^ i>r.rit \ 



3 

Since 0.5 of the mass of the KBO's is dust, we write that 



Mtotal 

Therefore, we find from equations (34) - (36) that: 

In 



form = 0.5^^^ (36) 



farm — 0.5 J- r- (37) 

In ) 

While the values of hmin and hmax are only poorly known, it seems for many comets that 
^min ~ 0.1 //m and hynax ~ 1 cm (see McDonnell et al. 1987, Manner et al. 1992). Results 
from meteorites (Kyte & Wasson 1986, Love & Brownlee 1995) and the Long Duration 
Exposure Facihty (Zolensky et al. 1995) also show that many particles larger than 100 /xm 
exist in the Solar system. With these extrapolations from the Solar System, we find from 
equation (37) that farm ~ 0.25. 



3.7. Infrared Excess 



We now elaborate upon the schematic discussion given in §2 to show that the infrared 
excess around a red giant, L excess i directly depends upon the rate of dust production from 
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the KBO's, Mkbo- Here, we include two additional physical processes not discussed in §2 
that control the grain's orbital evolution: stellar wind drag and grain-grain collisions. 

The stellar wind drag is straightforward to incorporate into the model. Re-writing the 
notation of Gustafson (1994) and Burns et al. (1979), the increase by "drag" in the inward 
drift velocity over that produced by the Poynting Robertson effect is given approximately 

by the factor (1 -|- MyjindC^ / L^) where My^ind is the stellar wind mass loss rate. In the Solar 
System, this "drag" created by the anisotropic recoil of solar wind particles off an orbiting 
dust grain typically scales as 0.3 of the Poynting Robertson drag (sec, for example, Gustafson 
1994). Below, we argue that the stellar wind drag around red giants is of similar magnitude. 

Consider now the orbital evolution of dust particles released from warmed-up KBO's. 
We assume that as seen from the host star, the dust ejected from the sublimating KBO's is 
confined to a solid angle, Q^bo- If the dust cloud is optically thin, the excess luminosity 
scales directly as the optical depth of the dust. If ndust,i denotes the density of dust particles 
of radius bi, we write that: 



M^* / '^bf ndust,i{D) dD (38) 



In equation (38), we assume that all the dust particles are formed at an initial outer distance 
from the star, Dinit, and they are all destroyed at some final, inner radius, D final- We also 
assume that the cross section of the grains is simply given by their projected surface area. 

We now show how the Poynting- Robertson and stellar wind drags affect the dust motion. 
Prom the equation of continuity, we write: 

where Ndust,i denotes the rate of production by number of dust grains from KBO sublimation 
and Vi{D) denotes the inward radial drift speed of the same grains. If Mdust,i denotes the 
rate of production by mass of the i'th type of particle, then: 

Ndusti = (40) 

Since the dust particles drift inwards under the action of the Poynting-Robertson drag, we 
write: 
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(Burns et al. 1979, Gustafson 1996). Note that since V — dD/dt and if we ignore the 
correction for stellar wind drag, equation (1) can be derived from equation (41). 

We are now able to compute the infrared excess from the system. Combining equations 
(38)-(41), we find: 

Prom §3.6, we write that 

. Mdust,i — MxEof orbit (43) 

Therefore, evaluating equation (42) and using equation (43): 

Lexcess = ^ MkbO f orbit ^1 + ^'"'^'^^ ^ In ^ ^^^J (44) 

Equation (44) is a fundamental result in this paper and shows that the excess infrared 
luminosity scales directly as the dust production rate, Mr Bo- 
lt should be recognized that equation (44) is a more complex version of equation (4) 
derived in §2. In equation (44), we employ Mrbo instead of Mkbo{^) /tRGB- Also, in 
equation (44) we include the wind contribution to limiting the particle's orbital lifetime in 
addition to the Poynting- Robertson effect. Finally, there are numerical constants that are 
different between equations (4) and (44). 

We now estimate the parameters required to use equation (44). First, we estimate the 
magnitude of the stellar wind drag. Currently, there are only a few first ascent red giants 
where the mass loss rate is measured with much confidence. Robinson, Carpenter & Brown 
(1998) report and M^ind for a Tau (K5 III) and 7 Dra (K5 III) of 394 Lq and -1.3 x 
lO"^""^ Mq yr"-*^ and 535 L0 and ~1.0 x 10^"*^^ M0 yr^-*^, respectively. These results yield 
values of Mwind(? / L^, of 0.49 and 0.28, respectively. For HR 6902 (02 II), the luminosity is 
550 Lq and Mwind is somewhere between 0.8 and 3.4 x 10^^^ Mq yr~^ (Kirsch, Baade & 
Reimers 2001). For this star, the correction to the Poynting Robertson drag is between 0.2 
and 0.9. Here, we simply assume that all red giants have the same stellar wind drag and 
adopt MwindC^/L* — 0.4, a value for this ratio close to that in the Solar System. 

To use equation (44), we also need to estimate D final /Dinn- In the absence of any other 
effects, the refractory dust would drift inwards until it reached a temperature of ~1000 K 
when it would be destroyed by sublimation. From equation (16) we can show that if T ~ 
170 K at Dinit, then Dfinai ~ 0.03 Dmit- However, before they sublimate, the particles' 
orbits might be disrupted as they drift inwards by grain-grain collisions. In such collisions. 
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the grains may shatter into pieces, some of which are small enough to be driven away from 
the star by radiation pressure. Shattering may also produce more smaller grains with a net 
increase in surface area of the particles. If the radial optical depth through the circumstellar 
dust is r, then, depending upon the velocity dispersion of the grains and their size and spatial 
distributions, the probability during one orbital period, P, of a grain-grain collision may also 
be of order r. Therefore, the mean time between collisions is P/r. With an orbital period 
at 45 AU of 300 yr and r = 10~^, the optical depth of interest here since in this case the 
infrared excess becomes detectable, the mean time between collisions of 3 x 10^ yr, which 
is comparable to the Poynting Robertson decay time for a dust particle of radius 100 /im 
orbiting around a red giant with a luminosity of 300 Lq. Therefore, as a very simple first 
approximation, grain-grain collisions are assumed to occur when the dust particle's orbital 
decay from the Poynting-Robertson drift has significantly changed its distance from the star. 
Although detailed calculations are required (sec, for example, Wyatt ct al. 1999, Krivov, 
Mann & Krivova 2000), we adopt D final = 0.5 Dinud^s a crude representation of the effects of 
grain-grain collisions. Our results for the infrared excess luminosity are not strongly sensitive 
to this value for Amt/-D/ma/ since, from equation (44), Ijexcess depends only logarithmically 
upon this ratio. 

Using the calculations for Mkbo described in §3.5 and equation (44), we show in Figure 
3 the results for Lexcess/L* as a function of for several models. We see that L^xcess/L* 
increases until it reaches its maximum value which occurs when Aa Models with 

initially some small KBO's and Ubreak — 1 km show a higher value of Mrbo than the model 
with only large KBO's and a^yreak — 200 km. Once the star's luminosity is high enough that 
Mkbo achieves saturation, then L^xcess is approximately constant and Lf.xcessj^'* decreases 
directly with L^,. From equations (30) and (44), the value for Lg^cess when Mkbo reaches 
saturation is: 




For the standard parameters adopted above, equation (45), an elaboration of equation (4) 
derived in §2, yields L^xcess ~ 0.02 Lq when T = 170 K. 

4. COMPARISON WITH OBSERVATIONS 

In §3, we have calculated the excess luminosity produced by a cloud of sublimating 
KBO's. We now show that this radiation may be detectable. 

To measure the amount of excess radiation that a red giant exhibits at 25 //m, we need 
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to determine its photospheric emission at this wavelength. Previous studies of first ascent 
red giants have shown that the bulk of the detected flux at 12 //m, F,^(12 /xm), emerges from 
a star's photosphere (see, for example, Jura, Webb & Kahane 2001, Knauer, Ivezic & Knapp 
2001). Here, we use the usual notation that v denotes frequency. We extrapolate from this 
measurement of the total flux at 12 /xm to estimate the photospheric contribution to F^(25 
/xm). From IRAS data. Jura (1999) found for the photospheres of the very brightest red 
giants that on average, F^(25 /xm)/Fj^(12 /xm) — 0.233. We define /ea:(25), the excess flux at 
25 /xm, by the expression: 

f (25) - ^ PA^^i^rn) _ 

/e.i^5j - ^ 233 ^^(12;,^) ^ ^^^^ 

In Figure 4, we show a histogram of /e2:(25) for 66 nearby red giants identified by a procedure 
described below. Over 90% of the red giants exhibit |/ea;(25)| < 0.10 implying that the IRAS 
data can be used to predict the photospheric flux at 25 /xm to within 10%. 

The histogram shown in Figure 4 is not symmetric around /ex (25) = 0. The sample of 
66 stars from which the data are obtained has a mixture of surface gravities and abundances, 
and the intrinsic distribution of /ea;(25) around its mean value may be skewed. The skewed 
histogram might also result if some of these 66 stars possess an excess at 25 /xm of ~ 5% 
of the photospheric flux. Better data should allow us to confirm the reahty of this hint of 
infrared excesses. 

The sample of stars is chosen from nearby red giants which are least subject to inter- 
stellar extinction. We select stars from the the Yale Bright Star Catalog with B-V > 1.0 
mag which also lie within 150 pc of the Sun (from Hipparcos data), have My < 0.0 mag 
and have 200 Lq < L < 300 < L0. We use 200 Lq as the lower bound of the luminosities 
that we consider because, as shown in Figure 3, for stellar luminosities lower than this vahic, 
the amount of the infrared excess becomes very sensitive to the initial numbers of KBO's 
smaller in radius than 200 km. We use 300 L0 as the upper bound for the stellar luminosity 
because for luminosities higher than this value, equation (44) begins to overestimate the 
infrared excess as Aa approaches a^ax- We compute the star's luminosity by assuming no 
reddening and applying the bolometric corrections from Flower (1996) for red giants with 
the assumption that MboI (Sun) = 4.74 mag (Bessell, CasteUi & Plez 1998). There are 66 
stars which satisfy these criteria and which were detected with IRAS. 

In our models described in §3, D final = 0.5 Dinu- Because this change in orbital radius 
is relatively small, we assume that the grains emit at a single temperature, determined from 
equation (16), by their initial distance from the star. Assuming that the infrared excess is 
characterized by a 170 K black body, then L^xcess ~ 1-5 vLi,{2b) where -L,^(25) denotes the 
specific luminosity at 25 /xm. 
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In Figure 3, in addition to showing the calculations for the standard model of L^xcessl L* 
vs. L^, we also plot the data for the 66 identified red giants. For each star, we use its distance 
and the non-color corrected flux from the IRAS Faint Source Catalog or, when that is not 
available, the flux from the Point Source Catalog, to estimate Li,{25). For the stars with 
/ex(25) < 0.1, we plot with a horizontal bar the quantity (0.1) (1.5)z/L;y(25 /im)/L* vs. L*. 
For stars with /ea;(25) > 0.10, we plot with a filled triangle the quantity /ea:(25)(1.5)i/Li,(25 
/im)/L^ vs. L^,. Of the six stars with /ea;(25) > 0.1, there are three stars, HR 736, HR 7956 
and HR 9029, which also have 60 //m excesses. The 60 //m excess for HR 9029 also has been 
noted by Plets et al. (1997). These infrared excesses are so small that it is very difficult to 
infer a dust temperature, and this evidence is at best a weak hint that these stars possess 
sublimating KBO's. 

From Figure 3, the we see that the data points usually he below the theoretical curve for 
the regime where initially there are many KBO's smaller than 200 km in radius. Therefore, 
we find that the production of dust from systems containing comet-like bodies can produce a 
detectable infrared excess if the total mass of the KBO's is at least 0.1 M^. We can also see 
from Figure 3 that if initially there are only KBO's with a radius larger than 200 km, then 
a KBO system similar to that in own Solar System would not produce a detectable infrared 
excess around most red giants. However, this calculation was performed with M^ig^O) = 
0.02 Mq. If, instead, we assume that Mbig{0) ?a 0.1 M®, then the predicted fluxes would 
be a factor of five larger than shown. In this case, the data points would then lie below 
the computed line and an infrared excess would be detectable. Thus, as anticipated in our 
discussion in §2, regardless of the exact size distribution, we expect that an infrared excess 
might be detectable if the total mass of the KBO system lying at 45 AU orbital radius is 
larger than 0.1 M®. 



5. DISCUSSION 

Our main argument is that as a star evolves on the red giant phase, it is possible to 
constrain the mass of the KBO's of the system. Here, we consider how our results compare 
with other investigations of the outer regions of our Solar System and its analogs among 
other stars. 

While sublimation of KBO's at an orbital radius of 45 AU might lead to a detectable 
infrared excess during a star's red giant evolution, our analysis does not provide strong 
constraints on the mass of KBO's much further than 50 AU. Melnick ct al. (2001), Ford 
& Neufeld (2001) and Ford et al. (2003) have proposed that in order to account for the 
detected circumstellar H2O and OH around the mass-losing carbon star IRC-(-10216, there 
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is a cloud of comets at ~300 AU with a mass of at least ~10 M®. In contrast, our typical 
upper limit to the mass in comet-like objects is ~0.1 M® at ~ 45 AU. This apparent mass 
gradient of KBO's is remarkable and deserves further study. 

Main sequence stars often possess from 0.01 M® to 0.1 M® of dust at orbital separations 
of typically between 10 AU and 100 AU (see Zuckerman 2001). Our results for evolving red 
giants are consistent with the hypothesis that there may be as much as an additional ~ 0.1 

in icy comets closer than ~50 AU around most main sequence stars. However, it appears 
that macroscopic KBO's with radii greater than 1 km do not overwhelmingly dominate the 
mass budgets of solid material orbiting beyond 10 AU around those main sequence stars 
with measured amounts of dust. 

The total mass of the KBO's in orbit around our Sun is somewhere between 0.01 M® 
and 0.1 Mq. Even the larger value is much less mass than predicted by extrapolations from 
the inner Solar System for the initial Solar Nebula (Luu & Jewitt 2002). Our upper limits 
for the mass of systems of KBO's around red giants shows that the analogs to the Kuiper 
Belt also have similarly "low" masses at orbital radii less than 50 AU. 

Very precise measurements of the fluxes from red giants in comparison with model 
atmospheres may allow for a more reliable determination of whether they possess infrared 
excesses. If the photospheric flux could be estimated to within 1%, instead of 10%, it would 
be possible to improve estimates of the masses of the KBO systems by a factor of 10. With 
such a capability, we could determine whether the hint in the data of some stars having 
KBO's of ~0.1 Me. 

In this paper, we have focused on the production of an infrared excess by sublimation of 
comets. Our calculations do not pertain to those red giants where an infrared excess with a 
characteristic dust temperature of ~70 K is measured that is produced by ~0.1 M® of dust 
at ~200 AU from the star (Jura 1999). The origin of this dust is uncertain; at least in some 
cases, confusion with interstellar cirrus is possible. Further work is required to understand 
the nature of this material. 



6. CONCLUSIONS 

We have computed the evolution of a system of KBO's when a star becomes a red giant. 
For KBO's at 45 AU from the host star, we flnd that the dust released during the sublimation 
of ice from these objects might produce a detectable infrared excess around red giants of 
luminosities in the range 200 Lq to 300 Lq if the mass of the KBO's is at least 0.1 Mq. To 
date, there is no strong evidence for such sublimating KBO's, but with improved data, it 
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may be possible to detect systems which are analogs to our Solar System's KBO's around 
first ascent red giants. 

This work has been partly supported by NASA. 
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Fig. 1. — We show a plot of theoretical calculations from Girardi et al. (2000) for the 
luminosity on the Red Giant Branch vs. time. The filled triangles show the results for a star 
of 1 Mq, Lq = 87 Lq and a starting time of 1.220 x 10^° yr while the open squares show the 
calculations for a star of 1.5 M0, Lq = 91 Lq and a starting time of 2.7742 x 10^ yr. The 
solid fine shows the prediction from equation (12) with tuoB — 14 Myr. The sohd curve hes 
within 20% of the detailed calculations for L < 1000 for both sets of models. 
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Fig. 2. — Plots of Mkbo from a cloud of KBO's vs. the star's luminosity. The solid line 
shows the calculations for the standard model consisting of both small and large KBO's with 
Dinit — 45 AU and Miargd^) = 0.1 while the dotted and short-dashed lines show the 
results for the cases where Anit — 40 AU and 50 AU, respectively. Note the approach to 
saturation for Mkbo given by equation (30). The long-dashed hue shows the calculations 
from equation (31) for a cloud of only big KBO's with Mug{Q) = 0.02 at Amt = 45 AU. 
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Fig. 3. — The solid curve shows a plot of Lexcess/L^, vs. for the standard model with 
both small and large KBO's described in the text with -Dj„jt = 45 AU and M;aj.ge(0) = 0.1 
M^. The dotted and dashed lines show calculations for the same parameters as the standard 
model except that we use Amt = 40 AU and 50 AU, respectively. The long-dashed line 
shows the results for the case with Dinu — 45 AU and only large KBO's with Mfeig(O) = 
0.02 Mq. We also display data for 66 stars identified by a procedure described in the text. 
The bars represent stars where /ex(25) defined by equation (46) is less than 0.1 while the 
triangles represent those stars where /ea;(25) > 0.1. 
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Fig. 4. — A histogram of the fractional excess at 25 /xm, /ea:(25), defined in equation (35), 
for the stars represented in Figure 3. 



